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[0001] This application is related to co-pending application serial number (Attorney 

Docket No. 1391-27200) titled "Acoustic Logging Apparatus And Method For Anisotropic Earth 
Formations." 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

[0002] Not applicable. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0003] The present invention relates generally to acoustic well logging. More particularly, the 
present invention relates to determining the velocity (slowness) of acoustic waves in an earth 
formation. 

Description of the Related Art 

[0004] In the art of acoustic logging, the formation properties of interest are one or more of the 
compressional wave velocity, the shear wave velocity and the Stonley wave velocity. These 
acoustic velocities are indicative of the formation's ability to hold and produce hydrocarbons. 
Operation of a typical acoustic formation tool may involve placing the tool in the well bore and 
firing one or more acoustic transmitters periodically, thus launching acoustic energy into the 
formation. The acoustic energy thus produced propagates along the formation wall in one or more 
propagation modes, e.g. compressional or shear wave modes. Receivers on the tool, spaced apart 
from the one or more transmitters, receive acoustic energy as the various waves move along the 
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wall past the receivers. The amplitudes and arrival times of the various acoustic signals at the 
receivers are indicative of the velocity of the acoustic waves (or slowness, being the inverse of the 
velocity) within the formation. 

[0005] Determining the acoustic velocity with early wireline acoustic logging tools involves an 
adaptation of data processing techniques used in seismic work. In particular, a method called 
semblance or coherence is used. United States Patent No. 4,594,691 to Kimball et al. (hereinafter 
the '691 patent) is exemplary of related art wireline acoustic logging tools that use this 
semblance/coherence measure for determining acoustic velocities within the formation. As 
exemplified in the '691 patent, determining the acoustic velocity using a coherence calculation is 
just a determination of how much two or more received waveforms resemble one another. The 
semblance/coherence determination itself, however, is not at all concerned with actual formation 
properties; rather, the power of the semblance/coherence measure is running the calculation on 
portions of each received waveform that should correspond, given the slowness of the formation. 
Since the formation slowness is the parameter of interest, the semblance/coherence measure is run 
multiple times at multiple slowness values, and the slowness values where the waveforms show the 
best semblance/coherence' are assumed to be the correct slowness values for the formation. 
[0006] Although the semblance/coherence calculation method has been relatively successful, the 
oil and gas industry as a whole has been searching for different, better and more efficient ways to 
calculate acoustic velocity in earth formations. United States Patent No. 5,541,890 to Tang 
(hereinafter the '890 patent) is exemplary of one such attempt. In particular, the '890 patent 
discloses that in the art of acoustic logging, any received signal may be synthesized by use of the 
other received signals. The '890 patent discloses synthesizing or estimating a received signal, and 
then comparing the synthesized signal to the actual received signal. For an array of possible 
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slowness values, the '890 patent discloses creating an object function which is indicative of the 
difference between the synthesized received signal and the actual received signal. In the '890 
patent, where that object function reaches a minimum (the difference between the synthesized 
signal and the actual signal are small), the values of the assumed slowness must match the actual 
formation slowness. Stated otherwise, the '890 patent records a series of waveforms at a plurality 
of receivers for each firing of the acoustic transmitter. The '890 patent then discloses synthesizing 
one received waveform from the remaining received waveforms, and comparing the synthesized 
received signal against the actual received signal. 

[0007] A similar technique is disclosed in U.S. Patent No. 5,081,611 to Hornby (hereinafter the 
'611 patent). The '61 1 patent discloses a method for determining slowness of an earth formation 
that involves calculating or estimating ray paths (ray tracing) and also estimating slowness values. 
The '611 patent discloses comparing the first arrivals (time of the first significant amplitude 
deflection in the received signal) predicted by the estimated ray path and slowness values to the 
actual first arrivals in each received signal. The patent discloses "back projecting" the error 
between the actual and calculated first arrivals (here agaio, in the received signals) and running the 
process of estimating ray paths and slowness values again until the error is small, thus revealing the 
actual ray paths and formation slowness. Much like the '890 patent described above, the '611 
patent is effectively estimating received signals, and comparing those estimated received signals to 
the actual received signals as part of determining slowness. 

[0008] While semblance/coherence may create visually pleasing results, determining slowness in 
this matter is not suitable for error estimation. Techniques described in the '61 1 patent and the '890 
patent disclosures may suffer from the same shortcomings. Consequently, an improved method to 
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determine the velocity or slowness of the various propagation modes in an earth formation is 
desired. 

SUMMARY OF SOME OF 
THE PREFERRED EMBODIMENTS 

[0009] The preferred embodiments of the present invention are directed to a method for 
determining the acoustic velocity of the various propagation modes of acoustic waves in an earth 
formation by assuming a transfer function of the earth formation and using the assumed transfer 
function to estimate a source signal based on each received signal. Preferably, the estimated source 
signals are compared to reveal an objective function from which the acoustic slowness of the earth 
formation may be determined. Preferably this is done for a plurality of assumed transfer functions, 
and minimas of the objective function indicate the formation parameters of interest. 
[0010] A first embodiment for comparing estimated source signals involves determining a variance 
of each of the estimated source signals against an average estimated source signal. This 
embodiment involves calculating an average of the estimated source signals, and then calculating a 
variance of the individual source signals against the average source signal. The variance values, 
which are the objective function in this case, may be plotted in such a way that the plot that looks 
similar to the time slowness plot used in semblance/coherence processing, but has two distinct 
advantages. First, the variance values have sharper transitions to the inflection points (minima in 
the preferred embodiments) than the semblance/coherence of the prior art, and therefore a more 
exact value of the slowness may be determined. Second, both of the variance values themselves, 
and the curvature of the objective function at the minima, are indicative of the accuracy of the 
slowness determination. That is, a particular plot of information may reveal a slowness value, and 
also indicates the relative accuracy of the slowness determination. 
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[0011] The second embodiment for determining the slowness of the formation of the preferred 
embodiments involves comparing each estimated source signal with its immediate neighbor to 
determine a difference, and then summing the square of the differences to obtain a value of the 
objective function. The values of the objective function may likewise be plotted in a time slowness 
plot to reveal a graph similar to a semblance/coherence plot familiar to one of ordinary skill in the 
art. The advantage of this second embodiment is that the plotted objective function has slower 
transitions from maxima to minima such that the determination of the minima points (determining 
the slowness values) may be accomplished with less computational overhead. In other words, 
determining inflection points in objective functions created by the differential comparison is at least 
as computationally fast as traditional semblance/coherence processing. Also, the values of the 
objective function at the minima and/or the curvature of the objective function at the minima is 
indicative of the error in the slowness determination. 

[0012] The disclosed device comprises a combination of features and advantages which enable it to 
overcome the deficiencies of the prior art devices. The various characteristics described above, as 
well as other features, will be readily apparent to those skilled in the art upon reading the following 
detailed description, and by referring to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0013] For a more detailed description of the preferred embodiments of the present invention, 
reference will now be made to the accompanying drawings, wherein: 
|0014] Figure 1 shows a wireline logging tool of the preferred embodiment; 
[0015] Figure 2 shows an exemplary set of received signals as well as an exemplary set of time 
slices used in the preferred embodiments; 
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[0016] Figure 3 shows an exemplary way to build a plot of the objective function in the preferred 
embodiments; 

[0017] Figure 4 shows an exemplary plot of a variance objective function at a single starting time; 
and 

[0018] Figure 5 shows an exemplary plot of a differential objective function at a single starting 
time. 

NOTATION AND NOMENCLATURE 
[0019] Certain terms are used throughout the following description and claims to refer to 
particular system components. This document does not intend to distinguish between 
components that differ in name but not function. 

[0020] In the following discussion and in the claims, the terms "including" and "comprising" are 
used in an open-ended fashion, and thus should be interpreted to mean "including, but not limited 
to...". Also, the term "couple" or "couples" is intended to mean either an indirect or direct 
electrical connection. Thus, if a first device couples to a second device, that connection may be 
through a direct electrical connection, or through an indirect electrical connection via other 
devices and connections. 

[0021] Although there may be some distinction in the oil and gas industry between logging- 
while-drilling (LWD) and measuring-while-drilling (MWD), throughout this specification and 
claims the terms are considered equivalent and mean generally any device or devices that makes 
measurements downhole (whether of formation properties, borehole properties, or properties of 
the drill string) during the drilling, logging and monitoring process. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENT 

[0022] Figure 1 shows an acoustic logging device 10 constructed in accordance with the preferred 
embodiments. In particular, Figure 1 shows a tool 10 disposed within a fluid filled borehole 12 at 
some distance below the surface 14. The tool 10 is preferably suspended within the borehole by 
means of a multi-conductor armored cable 16, and the cable 16 preferably couples to a surface 
computer (not shown). The tool 10 preferably comprises an acoustic transmitter 18 and a plurality 
of receivers 20A-D spaced apart from the transmitter 18 3 and also from each other. Although the 
preferred embodiment comprises four receivers, any number of receivers may be used without 
departing from the spirit and scope of the invention. 

[0023] Thus, the tool 10 of the preferred embodiment is a wireline device; however, the principles 
described herein may likewise be practiced in a measuring-while-drilling (MWD), logging-while- 
drilling (LWD) or any downhole measurement system. It must be understood that while the 
preferred embodiments of the present invention will be described in the context of a wireline 
device, the description in this manner should in no way be construed as limiting the general 
structures and processing methods disclosed herein to wireline devices. 

[0024] Operation of the wireline tool 10 involves placing the tool within a borehole 12 and 
allowing the tool 10 to drop or otherwise move below a formation of interest. The tool 10 is then 
slowly raised within the borehole, and while being raised, the transmitter 18 periodically fires, 
inducing acoustic energy into the formation. As the acoustic energy propagates through the 
formation, each of the receivers 20A-D receives a portion of the acoustic energy and converts the 
acoustic energy into electrical signals. The amplitudes of the received acoustic energy, as well as 
the time it took the energy to propagate along the distance between the transmitter 18 and the 
receivers 20A-D, are indicative of characteristics of the formation of interest. Determining these 
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characteristics involves analyzing the received signals to extract the information such as the 
acoustic velocity (or slowness, being the inverse of the velocity) of the various propagation modes 
such as the compressional wave, the shear wave or the Stonley wave. In the embodiment, where 
tool 10 is a wireline device, surface computers (not shown) are responsible for extracting this 
information. In an embodiment where the transmitter 18 and receivers 20A-D are located on a 
LWD or MWD device, making these acoustic measurements while the drill bit cuts through earth 
formations, some or all of the data processing required to extract the desired information may be 
done downhole. 

[0025] In broad terms, the processing method of the preferred embodiment comprises taking a 
plurality of received waveforms or received signals from tool 10, and calculating or estimating 
source signals or source wavelets based on each of the received signals by assuming a transfer 
function of the formation. The estimated source signal for each receiver or received signal is 
compared in some way with the other estimated source signals to create an objective function. For 
each set of received signals, this source estimation step is preferably performed multiple times, each 
time with a different assumed transfer function. Minimas of a graph of plotted objective function 
values are indicative of the slowness of the acoustic waves through the formation, and error of the 
calculation may be ascertained from the objective function graph as well. 

[0026] The source signal estimation will now be described in greater detail with reference to 
Figure 2. Figure 2 shows four received signals, 20A-D. Notice how the amplitude of the received 
signal of the closest receiver 20A is greater than the signal received by the more distant receivers, 
for example 20D, exemplifying attenuation of the acoustic signal as it propagates through the 
formation. Notice also how the waveforms shift out in time from the closest to the more distant 
receivers, exemplifying the finite speed of the acoustic waves within the formation. Determining 
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the slowness of the formation in the preferred embodiments involves assuming a transfer function 
for the formation. Using the assumed transfer function, an estimated source signal or source 
wavelet for each received signal is calculated. Estimating each source wavelet can be described 
mathematically as follows: 

W') = [7VT'tf,(<) (1) 
where S £STj is the estimated source signal calculated for the ith receiver, [TF] is the assumed 
transfer function of the formation for source to receiver propagation (and thus its inverse is used for 
estimating the source signal), R t is the ith receiver signal and t is time. For each received signal Rt, 
an estimate of the source signal S Esr . is created. 

[0027] As the above discussion implies, a primary component of the source signal estimation is the 
assumed transfer function [TF]. The assumed transfer function may be relatively simple, taking 
into account only the finite speed at which the acoustic signals propagate, or may be very complex, 
to include estimations of attenuation of the transmitted signal in the formation, frequency 
component dispersion, paths of travel of the acoustic signals, and the many different propagation 
modes within the formation (e.g. congressional waves, shear waves, Stonely waves). If desired, 
the transfer function used can even model the effects of the acoustic waves crossing boundaries 
between different layers of earth formations. For reasons of simplicity of the calculation, the 
preferred estimated transfer functions take into account only the propagation speed (slowness) of 
the acoustic energy in the formation. 

[00281 For a particular starting time, for example starting time Ti in Figure 2, and for a first 
assumed transfer function having an assumed slowness, portions of each received signal are 
identified as being related based on the transfer function. Stated another way, the portions of the 
received signals within the window or rectangular time slice 50 should correspond based on the 
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slowness of the assumed transfer function. The time width of the samples taken from each of the 
received signals (roughly the horizontal width of the rectangular time slice 50) is preferably at least 
as long as the source signal. In this way, an entire source waveform or source wavelet may be 
estimated. However, the time width of the samples taken from the received signals need not 
necessarily be this width, as shorter and longer times would be operational. While the preferred 
embodiment is to window portions of each signal as shown in Figure 3, and then estimate source 
signals for those portions of the waveforms within the window, windowing the received signals in 
this manner need not necessarily take place. It would be operable to estimate source signals using 
any portion of the received signals, and indeed each entire received waveform could be used to 
estimate a source signal, and this would still be within the contemplation of this invention. 
[0029] Thus, preferably, portions of the received signals in the rectangular time slice 50 are each 
used to create an estimated source signal. The estimated source signals are compared to create an 
objective function that is indicative of their similarity. The process of assuming a transfer function, 
estimating source signals or source wavelets based on received signals, and creating an objective 
function is repeated a plurality of times for each starting time. The rectangular time slices 50, 52 
and 54 are exemplary of multiple assumed transfer functions (each having a different assumed 
slowness) used in association with starting time Ti. Moreover, estimating source wavelets in this 
fashion (including multiple assumed transfer functions) is preferably repeated at a plurality of 
starting times within the received signals, represented in Figure 2 as Ti, T2 . . . Tn. 
[0030] The value of the objective function created for each assumed transfer function and starting 
time is preferably plotted in a graph as a function of the starting time and the slowness of assumed 
transfer function. As shown in Figure 3, the starting time of the source signal estimations, Ti, 
T2 . . . T n , are preferably the ordinate or X axis coordinate and the slowness is preferably the abscissa 
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or Y axis. For example, the value of the objective function calculated using portions of received 
signals 20A-D within rectangular time slice 50 is plotted at point 60 in Figure 4. Likewise, the 
value of the objective function calculated using portions of received signals 20A-D within 
rectangular time slice 52 is plotted at point 62, and the value of the objective function associated 
with the rectangular time slice 54 is shown at point 64. Thus, for each of a plurality of starting 
times, and for each starting time a plurality of slowness values associated with assumed transfer 
functions, a graph of the objective function is created. In the preferred embodiments, the value of 
the objective function plotted at each location is indicated by a color, with cooler colors (blues, 
purples) showing more similarity, and hotter colors (reds, oranges) showing less similarity, of the 
estimated source signals or source wavelets. The color schemes however are only exemplary and 
other schemes may be used without departing from the scope and spirit of the invention. . 
[0031] While comparing estimated source signals to obtain the objective function could take many 
forms, two embodiments will now be described ~ determining a variance of the estimated source 
signals, and summing a square of a difference of consecutive estimated source signals. 
[0032] Calculating the objective function of a first embodiment involves determining a variance of 
the estimated source signals. In broad terms, this embodiment involves calculating an average of 
the estimated source signals, and then calculating a variance using the average of the estimated 
source signals. In more mathematical terms, for each assumed transfer function, a series of 
estimated source waveforms or signals S KSTi are calculated using equation (1) above. From the 

estimated source signals, an average estimated source signal is preferably calculated as follows: 

5 ^0 = ^£w,(0 (2) 
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where is the average estimated source signal, N is the number of received signals, S ESTi is the 

source wavelet estimated for each received signal within the time slice, and t is time. 
[0033] The average estimated source signal is used to calculate a value representing the variance of 
the estimated source signals from the average estimated source signal. The variance is preferably 
calculated as follows: 

S 2 = Z(5 OT; (0-5 ra:t(?; (0) 2 (3) 
/=i 

where <> 2 is the variance. In this embodiment, the variance value is the objective function plotted in 
the slowness versus starting time graph (Figure 3). Large values of the variance indicate that the 
assumed transfer function (and its associated assumed slowness) did not significantly match the 
actual formation properties. Likewise, small values of the variance indicate that the assumed 
transfer function (and its associated assumed slowness) closely matched the actual formation 
properties. Thus, minimas of the objective function in the plot indicate the slowness of the 
formation, and the value of the objective function at those minimas is proportional to the error of 
the calculation. 

[0034] A second embodiment for calculating an objective function is based on determining a 
difference between each estimated source signal. In particular, and referring again to Figure 2, 
consider the portions of the received signals within rectangular time slice 50 associating with 
starting time Ti. As discussed above, using the assumed transfer function, an estimated source 
signal is created using the portions of the received signal within the time slice 50. In this 
embodiment, differences or differentials are calculated between each estimated source signal, for 
example between the source signal estimated from a portion of the 20A signal and the source signal 
estimated from the portion of the 20B signal. This difference is preferably calculated between each 
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succeeding estimated source signal, and the objective function in this embodiment is the sum of the 
square of each difference calculation. Much like plotting the variance objective function, this 
differential objective function is preferably plotted as a function of slowness and starting time, as 
shown in Figure 3. However, the graph obtained using the differential slowness calculation has 
slower transitions from maximas to minimas, which therefore makes determining the minima 
(indicative of the actual slowness of the formation) easier than in cases where the graph has 
relatively steep slopes between minima and maxima. More mathematically, the objective function 
of this second embodiment is calculated as follows: 

C= S(s EST+i -s ESTj ) 2 (4) 

i=i 

where C, is the objective function, and N is the number of received signals. Much like using the 
variance as the objective function, this differential objective function is preferably plotted in a 
slowness versus starting time graph with cooler colors representing less difference between 
received signals, and hotter colors representing greater differences. Known techniques may be used 
to determine minima of these graphs, and these minima are indicative of actual formation slowness. 
[0035] Referring now to Figures 4 and 5, there is shown exemplary plots of the two exemplary 
objective functions described in this specification. In particular, Figure 4 is an exemplary plot of 
the variance objective function against slowness. Figure 5 is an exemplary plot of the differential 
objective function against slowness. For purposes of illustration, consider that the variance value 
objective function plot of Figure 4 and the differential objective function plot of Figure 5 are 
based on the same set of estimated source waveforms. Each of the plots of Figures 4 and 5 could 
be, for example and without limitation, the values of an objective function for multiple assumed 
transfer functions along a single starting time. Referring briefly to Figure 3, the plots of Figures 
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4 and 5 could be the exemplary point 60, 62 through 64 associated with starting time Ti for the 
two methods of calculating the objective function of the preferred embodiments. 
[00361 E ac h of the objective function embodiments represented by the exemplary plots in 
Figures 4 and 5 5 have particular advantages. With reference briefly to the variance objective 
function of Figure 4 5 it has been described that the slowness of the formation may be determined 
by finding minimas of the objective function plot. Searching for minimas is typically done by 
computer programs. Computer programs, however, do not have the ability to see general trends 
in plotted information as does a human eye. In particular, a computer program searching for a 
minima in the variance objective function in the exemplary plot of Figure 4 cannot tell, 
comparing point 60 and 62, how close these points are to the minima of the plotted function. 
However, and referring to Figure 5, comparing points 60 and 62 of the differential objective 
function embodiment, a computer program can easily tell the direction and relative proximity of 
the minima. Thus, determining slowness in an embodiment using a differential objective 
function is faster than making the same determination in a variance objective function plot. 
Indeed, determining the minimas in a differential objective function plot rivals in computational 
speed the determination of standard semblance/coherence. Moreover, both the value of the 
objective function at the minima, for example point 66 in Figure 5, and the curvature of the plot 
of the objective function at point 66, are related to the error in the slowness determination. 
[0037] Computationally, determining the minima in the variance objective value system, such as 
that exemplified in Figure 4, is more difficult than in the differential objective function 
embodiment. However, the variance objective value embodiment has many advantages that 
outweigh any computational overhead required to find the minima (and therefore the slowness 
within the time slowness plot). In particular, and as exemplified in Figures 4 and 5, the variance 
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objective value embodiment has significantly sharper transitions from maxima to minima, which 
allow determinations of the formation slowness to be made more precisely than the differential 
objective function embodiment, and also more precisely than standard semblance/coherence. 
Further, both the value of the variance objective function at the minima as well as the curvature 
of the plot near that point are directly related to the error in the slowness determination. That is, 
a minima of a variance objective function values indicates by its location in the plot the slowness 
of the formation, but also the value itself shows the error. This error may also be determined as 
the second derivative of a function that describes the variance objective value near the minima. 
[0038] Numerous variations and modifications will become apparent to those skilled in the art 
once the above disclosure is fully appreciated. For example, the disclosed method for 
determining acoustic wave velocity and orientation may be implemented using any number of 
receiver levels and different receiver types for the acoustic logging tool. Throughout this 
discussion, the various earth formation characteristics were discussed with reference to finding 
minimas of the objective function; however, one of ordinary skill in the art could easily invert the 
values used, thus making a determination a search for maximum values in the plot, and this 
would be within the contemplation of the invention. Moreover, the two disclosed embodiments 
for finding the slowness need not be mutually exclusive — using both calculation methods in a 
single tool may be advantageous. Further, neither embodiment is preferred over the other. If the 
signal processing is taking place downhole, such as in a MWD or LWD system, the differential 
objective value method may be more advantageous. By contrast, if the tool of interest is a 
wireline tool, where computer speeds and memory size is of little concern, then the variance 
objective function may be more advantageous. It is noted, however, that either method may be 
implemented in either a wireline or MWD/LWD tool, and this would still be within the 
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contemplation of this invention. With regard to the differential objective function embodiment, 
calculating the objective function may equivalently be made by taking the absolute value of the 
difference, rather than the square. It is intended that the following claims be interpreted to 
embrace all such variations and modifications. 
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